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Abstract. Phase differences and amplitude ratios at different colour photometric bands are currently being used to discriminate 
pulsation modes in order to facilitate mode identification of Ar-driven non-radial pulsating stars. In addition to physical inputs 
(e.g., mass, T e ;f, etc.), these quantities depend on the non-adiabatic treatment of the atmosphere. This paper presents theoretical 
results concerning 8 Scuti pulsating stars. The envelope of each of these stellar structures possesses a convection zone whose 
development is determined by various factors. An interacting pulsation-atmosphere physical treatment is introduced which 
supplies two basic non-adiabatic physical quantities: the relative effective temperature variation and the phase lag (f> T , defined 
as the angle between effective temperature variations and radial displacement. These quantities can be used to derive the 
phase differences and amplitude ratios. Numerical values for these quantities depend critically on the a MLT parameter used to 
calculate the convection in the envelope. The dependence on a was analized and it was found that the use of colour observations 
may be of considerable importance in testing the MLT. Finally, examples are given of how a introduces uncertainties in the 
theoretical predictions regarding phases and amplitudes of photometric variations in 6 Scuti pulsating stars. 

Key words. Stars: 6 Scuti - Stars: atmosphere - Stars: oscillations 



1. Introduction 

Asteroseismology is presently being developed as an effi- 
cient instrument in the study of stellar interiors and evolution. 
Pulsational periods are the most important asteroseismologi- 
cal observational inputs. However, for non-solar like oscilla- 
tions, such as those characterizing 6 Scuti stars, the informa- 
tion contained in these periods is not sufficient for adequate 
constraints in theoretical predictions. Without additional obser- 
vational data, mode identification is, therefore, not feasible in 
non-solar like contexts. The primary obstacles are that the pul- 
sation modes of 6 Scuti stars: 1) do not lie within the asymp- 
totic regime and 2) may be affected by the "avoided crossing" 
phenomenon. Furthermore, rotation and eventual coupling also 
destroy any possible regular pattern. 

One way to obtain more information on the basis of photo- 
metric observations is to study the multicolor flux variations. 
The linear approximation to non-radial flux variations of a 
pulsating star was first derived by Dziembowski ( 1977b), and 
later reformulated by Balona and Stobie <1979al[T979b> and by 
Watson ( 1988 1. In Garrido et al. ( 1990l, the linear approxima- 
tion was successfully used to discriminate the spherical orders 
I of these pulsating star modes. Several attempts to fit observa- 
tions have shown that the method can be applied, at least to low 
rotational velocities (e.g., FG Vir in Breger et al. 1999 BI CMi 
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in Breger et al. l2TJfJ2l 4 CVn in Breger et al. [T9^9l VI 162 in 
Arentoft et al. 2001 and some other stars in Garrido 2000 for 
fast rotators, see Daszyhska et al. 2002). 

Most theoretical and numerical pulsation models have been 
developed using the adiabatic approximation (Christensen- 
Dalsgaard 119821 Tran & Leon 1995). However, pulsation is 
highly non-adiabatic in stellar surface layers, in which ther- 
mal relaxation time is either of the same order or even lower 
than the pulsation period. The accurate determination of the 
eigenfunctions in these layers will therefore require the use 
of a non-adiabatic description which includes the entire atmo- 
sphere. This procedure then makes it possible to relate multi- 
color photometric observables with such eigenfunctions. 

A number of authors (e.g., Dziembowski 1 1 977al Saio and 
Cox 1980 Pesnell 1990 Townsend 2002) have developed non- 
adiabatic codes and performed stellar pulsation calculations. 
However, these works have been carried out without a com- 
plete description of pulsation-atmosphere interaction. 

In this paper two approaches are compared. The first fol- 
lows Unno et al. (1989) for the equations and the numerical 
method used to solve them. Outer boundary conditions are im- 
posed at the level of the photosphere. This approximation will 
be referred to below as "without atmosphere". The second ap- 
plies a non-adiabatic pulsational treatment to the atmosphere 
derived by Dupret et al. (2002). This treatment has made it 
possible to obtain photometric observable values that are more 
realistic when compared with those generated by the "without" 
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atmosphere approximation (referred to below as "with atmo- 
sphere")- The latter approach allows theoretical predictions to 
be more directly connected with photometric observations than 
they previously had been by using only period comparisons. 

Theoretical and observed amplitude ratios and phase dif- 
ferences can be compared as observed in different color pho- 
tometric bands. Garrido et al. (1990 1 showed that the wave- 
length dependence of the limb darkening integrals is very weak 
for low I values. Thus, at least three color combinations, dis- 
tributed in the widest possible range of wavelengths, give con- 
sistent values for phase lag <p T and R. The latter parameter 
was defined by Watson (1988 1 and designed to measure de- 
partures from adiabatic conditions. Comparisons between non- 
adiabatic predictions and multicolor photometric observations 
have been made by Cugier et al. (1994i for (3 Cephei stars, 
Balona & Evers ( 1999) for 5 Scuti stars and Townsend (2002 1 
for Slowly Pulsating B stars. However, such comparisons were 
made without carrying out a detailed non-adiabatic pulsation 
treatment for the atmosphere. Photometric magnitude varia- 
tions have been reformulated by Dupret et al. ( 2003 ), thus mak- 
ing it possible to obtain theoretical non-adiabatic quantities 
which explicitly include the atmosphere (Dupret et al.|2002i. 
This reformulation has been employed in the study of [} Cephei 
and Slowly Pulsating B stars. 



In the present paper the approach of Dupret et al. ( 2002 
2003 1 is applied to 6 Scuti stars. Non-adiabatic theoretical mod- 
els "with" and "without" atmosphere, as well as their reso- 
lution algorithms, are presented. Solutions are compared in 
order to show the improvement introduced by including the 
pulsation-atmosphere interaction. One possible procedure for 
mode identifications using multicolor photometric diagrams is 
also suggested. For 6 Scuti stars, the non-adiabatic observables 
with which we are working can be directly related with the 
characteristics of their thin convective envelopes. These en- 
velopes are described by the Mixing Length Theory (MET) and 
parametrized by a, which is defined as the proportionality con- 
stant relating the mean path of a convective element with local 
pressure scale height. The a parameter can be constrained by 
searching for the best fit between theoretical and observed pho- 
tometric amplitude ratios and phase differences in color pho- 
tometric bands. Such observables also depend on the opacity 
variation in the Hell, Hel and HI ionization zones. 
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Table 1. Global characteristics of the equilibrium models 



out" atmosphere. These calculations then make it possible to 
derive (p T , |5T e ff/T e ff| and 6g e /g e , which are directly related 
with the photometric color variations. 

2.1. Non-adiabatic models "without" atmosphere 

Here the adiabatic and non-adiabatic equations have been de- 
rived following Unno et al. (I1989> . However, two variables 
regarding the stellar interior have been modified in order to 
directly link the interior and atmosphere pulsation equations. 
These variables are 
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where £ r is the radial displacement, P g , the gas pressure, cD, 
the gravitational potential, Lr, the radiative luminosity and 5X 
(resp. X'), the Lagrangian (resp. Eulerian) perturbation of the 
X variable. The other symbols bear their usual meaning. 

The linear non-adiabatic pulsation equations corresponding 
to these variables are described in Appendix A (Eqs. iA.H to 
( IA.6t ). and are obtained by neglecting the effects of rotation 
and the magnetic field. The "frozen" convection flux approxi- 
mation, 5Lc = and F' c = 0, has been chosen by following 
Unno et al. ( 1989 1. Boundary conditions are also taken from the 
same source. Internal conditions are described by Eqs. jA.9> to 
(A. 11) in Appendix A. Surface conditions are given here be- 
cause they are particularly relevant to the comparison of the 
"with atmosphere" (Sect. l2~2"t and "without atmosphere" treat- 
ment results: 



2. Theoretical models 

Equilibrium stellar models have been computed using the 
CESAM code (Morel Il997t . In accordance with this code, the 
stellar atmosphere can be described: 1) as a single layer (i.e., 
the photosphere) calculated in the Eddington approximation, or 
2) using the Kurucz equilibrium atmospheric models (Kurucz 
119931 to reconstruct the atmosphere from a specific Rosseland 
optical depth until reaching the last edge of the star. Table 
shows the global characteristics of the equilibrium models. In 
the convective core of these models, overshooting has been set 
to 0.2 times the local pressure scale height. 

The pulsational code begins by computing the adiabatic so- 
lution for fixed values of n and /, and this computation is used as 
trial input for the non-adiabatic computations "with" or "with- 
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Note that Eq. is obtained by perturbing the equation which 
defines the effective temperature (L R = AnR 2 cr rac i 77-), as well 
as by assuming that in the photosphere: 



ST = 5Teff 

T ~ T eS 



(5) 



2.2. Non-adiabatic models "with" atmosphere 

Eqs. ( IA. H - dA~6l are also used to describe the pulsation inside 
the star, whereas Eqs. jA.9> -(A. 11) describe the internal bound- 
ary conditions. However, the complete atmosphere pulsation 



equations are solved following Dupret et al. (2002i. The latter 
equations have been derived by taking into account that: 1) the 
radiation field is not isotropic in the atmosphere and the diffu- 
sion approximation is therefore not valid in this case; and 2) 
the radiation stress tensor cannot be represented by a diagonal 
matrix with a constant element: P R = (l/3)aT 4 . 

For these reasons, new approximations have been used with 
respect to those solved for the interior in order to obtain an im- 
proved description of the atmosphere. This new approach is 
based on the fact that the atmosperic thermal relaxation time 
is very short as compared to pulsation periods. Consequently, 
it can be assumed that during each pulsation period, the atmo- 
sphere remains in thermal equilibrium. Rewriting the equations 
given in Dupret et al. {2002 1, those corresponding to the atmo- 
sphere are: 
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and where a R is the radiative acceleration and r the Rosseland 
optical depth. The atmospheric equilibrium quantities for these 
equations have been obtained from the Kurucz model atmo- 
spheres (Kurucz ll993> . 

In comparison to those used for the "without" atmosphere 
treatment, the boundary conditions for these equations are: 

1 . The mechanical boundary condition is obtained by neglect- 
ing the derivative of the gas pressure perturbation at the 
stellar surface: 
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3. The thermal boundary condition, following Dupret et al. 
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A connecting layer is required which marks the border be- 
tween the interior (Eqs. IA.H to JA.6I >) and atmosphere treat- 
ments (Eqs. (jfji to Jilt ). The boundary conditions for the equi- 
librium models are calculated at the same connecting layer 
by imposing an adequate match with the Kurucz model atmo- 
sphere (Kurucz fl993> . Two restrictions are imposed in order 
to choose the location of this layer: 1) the diffusion approxi- 
mation must be valid and 2) the layer must be located outside 
the convective envelope, since the atmospheric non-adiabatic 
treatment is not justified inside a convection zone, and because 
consistency requirements exist between the interior and atmo- 
sphere treatments. 

One of the uncertainties of the "with" atmosphere global 
algorithm is the location of the optical depth in which this 
connecting layer is defined, being that, despite the above- 
mentioned constraints, a small set of possible locations still ex- 
ists. Fig^shows the non-adiabatic results for two 2 Mq stellar 
models at different evolutionary stages, locating the connect- 
ing layer at r = 1, 1.1, 1.2 and 1.3. Note that these results are 
not significantly affected by this choice. The connecting layer 
is selected at r = 1 . 

The two additional variables ST e ff and 5g e can be derived 
from: 
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These quantities are determined at the photosphere (log T = 
log T e ff ). The results are not significantly sensitive to the choice 
of this layer: a slight change of location implies changes which 
are always smaller than 5 percent for |<5T e ff/T e ff|, smaller than 1 
degree for the phase-lag <p T and which are negligible for 6g e /g e . 

3. Multicolor photometry 

An important application of our non-adiabatic code is that its 
theoretical predictions can be tested against multicolor photo- 
metric observations. The non-adiabatic quantities |(5T e ft7T e ff|, 
5g e /g e and 4> T can be related to the photometric observables: 
amplitude ratios and phase differences between filters. In a 
one-layer linear approximation, the monochromatic magnitude 
variation of a non-radially pulsating star is given by: 
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Fig. 1. Non-adiabatic quantities |<5T e ff/T e ff| (top) and <p T (bot- 
tom), as a function of the pulsation constant Q (in days) for dif- 
ferent modes with spherical degrees I = 0, 1, 2, 3. Two 2.QMq 
models are solved at different evolutionary phases, with a MLT 
parameter a — 1 and the CEFF equation of state. Different val- 
ues correspond to different locations of the connecting layer: 
t = 1, 1.1, 1.2, 1.3. "+" are for the model with X c = 0.25 and 
"x"forX e = 0.55. 
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fj and f g are the relative amplitudes of local effective tempera- 
ture and gravity variations for a normalized radial displacement 
at the photosphere, <p T is the phase difference between the rel- 
ative effective temperature variation and the relative radial dis- 
placement and cr is the pulsation frequency, a corresponds to 
the real amplitude of the relative radial displacement, P'" is the 
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associated Legendre function and i the inclination angle be- 
tween the stellar rotation axis and the observer line of sight. 

One appropriate way to test multicolor theoretical predic- 
tions is to construct phase-amplitude diagrams corresponding 
to well-chosen combinations of photometric bands. In these di- 
agrams the theoretical results corresponding to modes of dif- 
ferent spherical degrees I occupy well differentiated regions. 
This enables the identification of t by searching for the best 
fit between theory and observations. As is shown in Sect. |4] 
the non-adiabatic quantities /j, f g and cp T , which play a major 
role in Eq. ( I19> . are highly sensitive to the characteristics of the 
convective envelope. The confrontation between the theoretical 
and observed amplitude ratios and phase differences can thus 
be used to constrain the physical conditions of this convection 
zone. 



4. Applications to 5 Scuti stars 

As a prerequisite to the presentation of the "with" atmosphere 
treatment results, of particular interest is that the differences be 
analized between the non-adiabatic observables appearing as a 
consequence of "with" and "without" atmosphere descriptions. 
In Fig. |2] these differences are displayed for a 1.8Mq model 
in the middle of the evolutionary phase (X c = 0.44) towards 
the exhaustion hydrogen in the core. For the relative ampli- 
tude of the local effective temperature variation, the effect of 
introducing the atmospheric treatment is slight, though signifi- 
cant, especially around the overstable modes (between Q=0.02 
and Q=0.012 days in this model). However, the phase lag <p T 
is more sensitive to this treatment for all frequencies and dis- 
plays a mean difference of 30° for the overstable modes. As of 
the latter modes, these quantities increase with the frequency. It 
should be noted that the differences in 1 — 1 between the "with" 
and "without" atmosphere models are negligible. Furthermore, 
it should be stressed that the frequencies obtained using both 
treatments are almost identical, i.e., the only difference be- 
tween the two treatments appears in the outermost layers, and 
this has no effect on the frequency results. In Fig- E] it can be 
seen that the non-adiabatic results for a fixed model are inde- 
pendent of the spherical degree / forp-modes. 

Only pulsation models "with" atmosphere will be consid- 
ered in what follows. In Fig. non-adiabatic results are com- 
pared as obtained for models of 1 .8Mq and X c = 0.44, and with 
different values for the MLT parameter (a = 0.5, 1, 1.5). The 
non-adiabatic quantities | <5T e ff /T e tf| and <p T are significantly af- 
fected by the value of a, especially the phase lag. A more de- 
tailed description of the sensitivity of the non-adiabatic results 
to a is given in Fig. @] This figure displays the ratio between 
radiative and total luminosity (top), the convective efficiency 
(middle) and the luminosity phase lag (pi = 0( ^) ~ 4\r) ^ ot " 
torn panel), all three as a function of the logarithm of tem- 
perature. Note that <p L closely follows <p T at the photosphere 
throughout Eq. d!7i . 

The bottom panel of Fig. |4] shows that along the radius of 
the star there are two zones in which a phase lag is introduced: 
the first located in the partial ionization zone of Hell and the 
second, in the surface convection zone (partial HI and Hel ion- 



1.8 M , oi=l, CEFF, X c =0.44 




0.01 0.015 0.02 0.025 

Q (in days) 

1.8 M , ct=l, CEFF, X c =0.44 



11.03 



0.035 




0.01 



0.015 



0.02 0.025 
Q (in days) 



0.03 



0.035 



Fig. 2. Non-adiabatic quantities |{5T e fF/T e ff| (top) and <p T (bot- 
tom), as a function of the pulsation constant Q (in days) for 
different modes with spherical degrees I = 0, 1,2,3. A model 
of a I.8M0 is solved with X c = 0.44, a MLT parameter a - I 
and the CEFF equation of state. Results obtained "with" (+) 
and "without" atmosphere (x) in the non-adiabatic treatment 
are compared. 

ization). In this panel it can also be seen that the the phase lag 
sensitivity to a appears in the surface convection zone. 
The convective efficiency in MLT, defined as 
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(Cox 119801 1. and the convection zone temperature gradient 
are both directly linked to a. Phase-lags originate in the en- 
ergy conservation equation due to the introduction of explicitly 
imaginary parts. For a radial mode, and freezing the convective 
luminosity, the equations are: 
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where M r is the mass enclosed in a sphere of radius r, and 
8L R 6k J r A 6T d^f/dlnr 
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Fig. 3. Non-adiabatic quantities [5T e ff/T e ff| (top) and <\> T (bot- 
tom), as a function of the pulsation constant Q (in days) for 
different modes with spherical degrees I — 0, 1, 2, 3. Models of 
I.8M0 are solved with X c = 0.44 and with different values for 
the MLT parameter: a = 0.5, 1, 1.5. Note the sensitivity of the 
non-adiabatic results with respect to a. 



Phase-lags originate mainly through the interplay between 
the different terms of Eq. ( I26> , which affect the right hand side 
of Eq. ( I25t only when the thermal relaxation time is sufficiently 
small. A first phase-lag is introduced in the partial ionization 
zone of Hell (bottom panel of Fig. 0). The partial ionization 
produces an opacity bump and a considerable decrease in the 
adiabatic exponents. This significantly affects 6k/k. A second 
source of phase lag occurs in the surface convective zone (par- 
tial ionization zone of HI and Hel). In this case variations in 
a primarily influence the temperature gradient and the size of 
the convective zone. The latter two then affect the phase-lag 
throughout the above equations. 

As is shown in Fig.^ the non-adiabatic results are highly 
sensitive to the evolutionary phase of the star. As the star 
evolves in the HR diagram, changes in these quantities are 
larger for the relative variations of the effective temperature 
than for the phase lags. Changes in the evolution phase do not 
qualitatively affect the non-adiabatic results in the same manner 
as do changes in a. In the latter case, only the characteristics of 
the convective zone are influenced. In contrast, changes in the 




1.8 M () , CEFF, X c =0.44 




Fig. 4. Radiative luminosity over total luminosity (top panel), 
convective efficiency (middle panel) and the luminosity phase 
lag (p L = ^(x) _ ^(t") (bottom panel) as a function of the loga- 
rithm of temperature, for the fundamental radial mode. 1.8Mq 
models are solved with X c = 0.44, and three different values 
of the MLT parameter: a = 0.5, 1, 1.5. Note the differences ap- 
pearing in the superficial convection zone as a values rise. 



evolution phase have an effect on both the characteristics of the 
superficial convective zone and the position of the Hell partial 
ionization zone. 

Fig.|5]is obtained by plotting the phase lag and |<5T e ff/T e ffl 
versus logT e ff for the fundamental radial mode of two com- 
plete tracks of 2.0M o and 1 .8M stars. Three values of the MLT 
parameter a (1.5,1.0 and 0.5) have been used. When stellar 
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in two complete tracks of 2.0M o and 1.8M Q stars for three dif- 
ferent values of the MLT parameter a =1.5,1 and 0.5 
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Fig. 6. /c-phase <p K (top) and convective-phase <f> c <mv (bottom) as 
a function of log T e fj for the fundamental radial mode in two 
complete tracks of 2.OM and 1.8M stars for three different 
values of the MLT parameter a =1.5,1.0 and 0.5 



temperatures are high enough to originate a negligible external 
convective zone, the values are independent of a\ this only oc- 
curs for the 2.0M Q models presented. For cooler models, how- 
ever, convection becomes more efficient, thus producing differ- 
ent phase lags and relative effective temperature variations for 
similar effective temperatures (Balona & Evers 1999 1. 

Fig. prefers to the behaviour of these phase lags. Each step 
in the figure makes it possible to view the phase lag, <p T , as 
the sum of two phases (0l = </>k + </>com>)- The latter can be 
referred to as the K-phase (<p K = phase difference between 180° 
and <pi (at log T = 4.5)) and the convective-phase (<p com , = 
phase difference between (pi (at log T = 4.5) and <Pl (at log T = 
log T e ffJ). Fig.|6]is obtained by plotting both versus log T e ff- 

The A:-phase lag <p K displays a behaviour which is directly 
related to the position of the model in the HR diagram and 
which is independent of the a parameter. Therefore, <p K shows 
the contribution of only one of the phase sources, the opacity 
variation, since convection is not efficient in this part of the star. 

The contribution to the convection zone phase, in this case 
4>com; is in turn a sum of two phase sources. For hot models 
the k driving mechanism is not yet efficient enough to intro- 
duce a phase lag, and the behaviour of <f> conv depends mainly 



on the convection treatment. When the convective layers are 
not efficient, models with different values of a produce the 
same values of <p conv . At log T e ff * 3.91 the convective lay- 
ers become sufficiently efficient to distinguish between models 
using different values of a, independently of the mass. From 
log T e ff - 3.88 to cooler models, an increase is observed in the 
value of this phase lag, again displaying different behaviours 
for different masses and values of a. This is accounted for by 
the fact that, within this range of temperature, the k driving 
mechanism in the HI and Hel ionization zone becomes efficient 
enough to make a significant contribution to <p conv . 

Fig.^is obtained by displaying |(5T e ff/T e jf| versus the inte- 
gral of the convective efficiency, as defined in eq. ( 1241 for all 
models. The behaviour of these quantities can be distinguished 
independently of the mass of the model, which in this case is 
exclusively a function of the evolution phase and a. 

5. Applications for multicolor photometry 

All of the above-mentioned calculations have a direct influence 
on the phase difference - amplitude ratio diagrams. Now phase 
lags, as well as relative variations in |<5T e ff/T e fr| and 5g e /g e 
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can be used to overcome the uncertainties in previous phase- 
ratio color diagrams. In Garrido (2000) these discrimination 
diagrams were made using parametrized values for departures 
from adiabaticity and phase lags. The only remaining degree 
of freedom is now the choice of the MLT a parameter in order 
to describe the convection. Therefore, discrimination diagrams 
depend only on this parameter, as is shown in Fig.[8]where the- 
oretical predictions are plotted for two specific Stromgren pho- 
tometric bands ((b — y) and y) with a given theoretical model 
using three different MLT a parameters in the fundamental ra- 
dial mode regime (pulsation constant near 0.033 days) and in 
the 3rd overtone regime (near 0.017 days). 

A clear separation between the Z-values exists for periods 
around the fundamental radial. Similar behaviour is found for 
other modes in the proximity of the 3rd radial overtone, al- 
though for these shorter periods some overlappings start to ap- 
pear at the lowest /-values. They also show the same trend as 
in the fundamental radial mode: high amplitude ratios for low 
MLT a and the spherical harmonic I — 3. A more detailed de- 
scription of these diagrams using different models and photo- 
metric bands will be given in a forthcoming paper. 
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Fig. 8. The top panel shows theoretical predictions for two spe- 
cific Stromgren photometric bands ( (b - y) and y) for a given 
theoretical model using three MLT a parameters in the fun- 
damental radial mode regime (pulsation constant near 0.033 
days). The 3rd overtone regime (pulsation constant near 0.017 
days) is shown in the bottom panel. 



6. Conclusions 

This paper uses the CESAM package to generate equilibrium 
models. A new non-adiabatic pulsation code has been devel- 
oped and applied to the study of 6 Scuti stars. This code 
takes into account the stellar atmosphere in two ways: 1) as 
the boundary layer of the star ("without" atmosphere); in this 
case photospheric observables such as amplitudes and phases 
of both T e ff and gravity, reflect the boundary conditions we 
impose; and 2) particular consideration is given to pulsation 
treatment ("with" atmosphere), thus extending the star beyond 
the photosphere. As a consequence of applying the "with" at- 
mosphere treatment, photospheric observables become deter- 
minable as a solution to a set of differential equations. These 
equations have been described by Dupret et al. (2002 1. 



The non-adiabatic code presented in this paper enables 
the determination of the photometric observables \6T e ff/T e ff\, 
\5g e lg e \ and <p T . The results of the "without" atmosphere ap- 
proach strongly depend on the choice of external boundary con- 
ditions. In contrast, the "with" atmosphere approach only im- 
plies those uncertainties brought about by the physical assump- 
tions imposed in order to obtain the differential equations. 

The comparison of the results generated by both treatments 
shows that the values for \ST e ff/T e ff\ are relatively similar 
in most of the spectrum. However, differences become no- 
table for shorter periods, in which "without" atmosphere solu- 
tions are larger than those of the "with" atmosphere treatment. 
Differences in <p T are significant and have a mean value of ap- 
proximately 30° within the range of the periods studied. The 
non-adiabatic results presented here are highly sensitive to the 
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characteristics of the superficial convective zone, parametrized 
by using the mixing length parameter a. In particular, we have 
shown that there are two regions in which the phase lag orig- 
inates. A first phase lag takes place in the partial ionization 
zone of Hell, where the k mechanism drives the oscillations. 
This phase lag is very sensitive to the evolution phase. A sec- 
ond one occurs in the convective envelope (partial ionization 
zone of HI and Hel) and is sensitive mainly to a. Though to a 
lesser extent, this phase lag is also sensitive to the evolutionary 
phase because of the k driving mechanism and the size of the 
convection zone change. 

In this paper it is shown that theoretical photometric am- 
plitude ratios and phase differences are very sensitive to non- 
adiabatic calculations. The treatment and location of the outer 
boundary conditions, here referred to as "with" and "with- 
out" atmosphere, have been shown to give significantly dif- 
ferent theoretical predictions. Consequently, our improved 
non-adiabatic treatment of the atmosphere will probably en- 
able more accurate photometric mode identifications. Our 
results confirm the recent theoretical results of Daszynska- 



Daszkiewicz et al. ( 2003 1 which show that phases and am- 
plitudes in different colours are affected by the choice of the 
a parameter, although the authors calculated the model atmo- 
spheres in the Eddington approximation. By determining the 
best fit between theory and observations, it may thus be possi- 
ble to constrain the MLT parameter a. 
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Appendix A: Stellar interior equations and inner 
boundary conditions 

These equations have been developed following Unno et al. 
( 1989 1. Eigenfunctions + and ^- have been replaced 

by -s^ and ^ to adapt the equations to the atmospheric ones. 
This was done to make the transition in the connecting layer as 
smooth as possible. Having introduced these modifications, the 
nonadiabatic nondimensional equations for the stellar interior 
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where jc = r//?, and R is the photometric radius. The rest of the 
variables follow the definitions given in Unno et al. {1989 1: 
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N is the Brunt-Vaisala frequency and, M the total mass of the 
star, k, is the opacity and, sn> is the energy generation rate. 

The interior boundary conditions are common to both pul- 
sational treatments, being that the stellar interior is solved using 
the same equations. In our formalism, these conditions are: 

1. The mechanical boundary condition 
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3. The thermodynamical boundary condition 
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